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ABSTRACT 

We present and interpret new spectropolarimetric observations of the mag- 
netic white dwarf WD 1953-011. Circular polarization and intensity spectra of 
the Ha spectral line demonstrate the presence of two-component magnetic field in 
the photosphere of this star. The geometry consists of a weak, large scale compo- 
nent, and a strong, localized component. Analyzing the rotationally modulated 
low-field component, we establish a rotation period Prot = 1.4480 ± 0.0001 days. 
Modeling the measured magnetic observables, we find that the low-field compo- 
nent can be described by the superposition of a dipole and quadrupole. According 
to the best-fit model, the inclination of the stellar rotation axis with respect to 
the line of sight is i ~ 20°, and the angle between the rotation axis and the 
dipolar axis is /3 ~ 10°. The dipole strength at the pole is about 180 kC, and the 
quadrupolar strength is about 230 kC. These data suggest a fossil origin of the 
low-field component. In contrast, the strong-field component exhibits a peculiar, 
localized structure ("magnetic spot") that confirms the conclusions of Maxted 
and co-workers. The mean field modulus of the spot {\Bspot\ = 520 ± 7 kC) to- 
gether with its variable longitudinal magnetic field having a maximum of about 
-1-400 kC make it difficult to describe it naturally as a high-order component of 
the star's global poloidal field. Instead, we suggest that the observed strong-field 
region has a geometry similar to a magnetic flux tube. 
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1. Introduction 



At present, there are more than one hundred known isolated magnetic white dwarfs 
(MWDs) wi th magnetic field 
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2007 
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2003 
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2004 
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20071 ). It is generally assumed that the magnetic fields of the strong-magnetic MWDs (those 
with MG-str ength fields) ar e organized as low-order multipolar fields with dominating dipolar 
components (|Putneylll999l ). The rotation periods and surface mag netic fields of the strong- 
magn etic MWDs are believed to be stable on long time scales (jSchmidt fc Northworthy 
19911 ). suggesting that their fields are fossil remnants of the fields of their progenitor stars. 



A comparison of the field strengths and incidence statistics of the strong-mag netic MWDs 



with magnetic fields of Ap/Bp main sequence stars support this assumption (lAngel et al. 
198lh . 



Despite the progress with the strong-magnetic MWDs, the magnetic properties of the 
weak-field degenerates are only poorly know n. Presently, only a few white dwarfs with 
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20071 ). Their rota tion and field g eomet ries a re poorly studied, although some 



progress has been achie ved by 



Maxted et all (120001) and IWade et all (120031) with a study 



of WD 1953-011 and bv IValvavin et al.l (120051 ) with WD 0009+501. iMaxted et al.l (I2OOOI ) 
established that the magnetic morphology of WD 1953-011 can be described by both low-field 
{B ~ 90 kG) and strong-field {B ~ 500 kG) components. Some evidence for the presence of 



a non - dipolar (quadrupolar) component was also found in WD 0009+501 by IValyavin et al. 



Motivated by the results of lMaxted et al.l (|2000l ). we have undertaken collaborative spec- 
tropolarimetric monitoring of WD1953-011. In this p aper we report resu lts o f these observa- 



tions and analyse them in the manner presented by I Wade et al.l (|2003l ) and IValyavin et al. 



( I2OO5I ). Our goal is to determine precisely the magnetic morphology of this degenerate. 
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2. A few preliminary remarks 



O ur preliminary an alysis of the spectropolarimetric data obtained with FORSl at the 
VLT (I Wade et al.l l2003l ) revealed significant variability of the Stokes / and V spectra of 
WD 1953-011 due to rotation, with a period of about 1.45 days. In Fig. [1] we show those 
results which illustrate the variation of the Stokes /, V, Q and U profiles of the Ha line with 
rotational phase (phase increases from top to bottom in the figure). As can be seen, the 
Stokes / profile is strongly variable. The central S-wave of the Stokes V profile is almost 
constant. However, near rotational phase 0.6 additional broadened Stokes V signatures ap- 
pe ar in the Ha w i ngs. T hese signatures correspond to the weak "satellite features" observed 
by iMaxted et al.l (120001 ) in the wings of the Ha profile at these phases, and were attributed 
to the presence of a high-field magnetic structure. Linear polarization Stokes Q and U 
signatures are only marginally detected at several rotational phases. 

In this pap er we extend this a nalysis using addition observational material obtained 
with the AAT (IMaxted et al.l l2000l ) and with the 6-m Russian telescope BTA, and using 
more sophisticated modeling techniques. 



Observations 



Spectropolarimetric observations of WD 1953—011 were obtained in service mode be- 
tween May and June 2001 with FORSl on the ESO VLT. FORSl is a multi-mode instru- 
ment for imagi ng and multi-objec t spec troscopy equipped with polarimetric optics, and 
is described by lAppenzeller et al.l (119981 ). For this work, FORSl was used to measure 
Stokes IQUV profiles of WD 1953—011 at 12 different rotation phases, using grism 600 R 
(plus order separation filter GG435), which covers the interval 5250 A- 7450 A. With a slit 
width of 0.7", the spectral resolving power was about 1650. To perform circular polarization 
measurements, a A/4 retarder wavep late and a WoUaston prism are inserted in the FORSl 
optical path (see lAppenzellerl (119671 )). The A/4 retarder waveplate can be rotated in 45° 
steps. 

To perform linear polarization measurements, a A/2 retarder waveplate is used, which 
can be rotated in 22.5° steps. At each epoch. Stokes V was measured taking two 420 s 
exposures: one with the A/4 retarder waveplate at —45°, and one with the A/4 retarder 
waveplate at +45°. Stokes Q and U were measured taking four 600 s exposures with the 
A/2 retarder waveplate at 0°, 22.5°, 45°, and 67.5°. The Stokes V/I circular polarization 
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spectrum was then obtained by calculating 

V 1 



V 



(r_45 - r+45) 



where 



r-r 



(1) 



(2) 



In Eq.(2) /° is the flux measured in the ordinary beam and is the flux measured in the 
extra-ordinary beam, obtained with the A/4 retarder waveplate at angle a. Similarly, the 
linear polarization was obtained by calculating 

Pq = 1 = |(^o-r45) 

Pu = J = 2 (^"22. 5 — ^67.5) 

where is defined by Eq. [2], except that (3 refers to the position angle of the A/2 retarder 
waveplate. Fluxes /° and were obtained from the raw data after bias correction and 
wavelength calibrations performed using standard IRAF routines. 

These observations were supported by a short observing run at the 6-m Russian telescope 
BTA where we obtained additional J, V series of spectra of WD 1953-011 using the UAGS 
spectropolarimeter, with nearly the same re solving power as i n the observati ons with FORSl. 



he i nstrument is described in detail by lAfanasief et al.l (119951 ) and by iNaydenov et al. 



( 2002). The observa t ional techniqu e and data reduction are similar to those described by 



Bagnulo et al.l (120001 . |2002| ) and by IValyavin et al.l (120051 ) . A comparative analysis of the 
spectropolarimetric data obtained from the different telescopes showed identical results. A 
comparison of the Stokes V spectra obtained with the VLT and BTA is illustrated in Fig. [2J 
The spectra are obtained at different times but similar rotational phases. 

In addition to the spectropolarimetric data from the VLT and BTA, in this paper we 
also use high-resolut ion spectroscopic data (Stokes /) obtained at the AAT and described by 



Maxted et al. 



J2000h. Together with the spectropolarimetry, these data extend the analysis 



presented by lWade et al.l (120031 ) to a much longer time base. Table [T] gives an overview of all 
the observations. In the table: JD is the Julian Date; Exp is an equivalent exposure time 
of an observation; Stokes is the observed Stokes parameter {I,V,Q or f/); and Telescope 
is telescope used (VLT, AAT or BTA). 



4. Mean field modulus of the magnetic field of WD 1953-011 

We begin with an analysis of the low and high-resolution Stokes I spectra, extending 
over 5 years. These spectra are used for establishing the rotation period of the star, as well 
as the mean field modulus of the low- and strong-field components. 
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4.1. The low-field component 



The low- field component of WD 1953-011 wa s first discovered spectroscopically by lKoester et al 

( 1998 ) and described in detail by Maxted et al. ( 2000 ). The mean field modulus, |-Bg0 



ex- 



hibits a low-ampl itude variation due t o the star's rotation, with a period estimated between 



hours and days (IMaxted et al.l |2000| ) . These conclusions were made on the basis of the 



high-resolution spectroscopy of the Zeeman pattern in the Ha core. 

In our low-resolution FORSl and BTA observations, Zeeman splitting attributed to the 
low-field component cannot be resolved spectroscopically. In the spectra the splitting is re- 
vealed as an additional variable broadening and desaturation of the Ha core. In this case, 
measurements of the low-field component can be carried out by an analysis of the equivalent 



width s (EWcore) of the Ha core. Using field modulus l-B^I values determined by lMaxted et al. 



( 120001 ) from an analysis of individual high-resolution Ha line profiles and measuring equiva- 
lent widths of the Ha cores, we may try to calibrate the relationship EWcore^ l^cl to allow 
us to determine \Bg\ in the low-resolution spectra. To obtain the required calibration, we 
estimated equ i valen t widths of the Ha core from the high resolution spectra obtained by 



Maxted et al.l (120001 ). In order to work with measurements having a uniform resolution, all 
high-resolution spectra were convolved with a gaussian instrumental profile to reproduce the 
spectral resolution of FORSl and UAGS. The resultant spectra are presented in Fig. [31 

As one can see in Fig. [3l the profiles are strongly variable. The central intensity of the 
core also correlates with the intensity of the strong-field Zeeman features which are found in 
the wings of the Ha profile (Fig. [3j the two satellite features at ±10 A around the Ha core). 
This correlation (the higher the intensity of the features, the weaker the central intensity) is 
due to the fact that the spot, which appears periodically on the visible disc due to rotation, 
redistributes the fiux according its projected area. It is also seen that the width of the Ha 
core is variable itself due to the variable Zeeman pattern of the line core. Therefore, to 
minimize the infiuence of the variable high-field spectral features in measurements of the 
central Zeeman pattern attributed to the low- field component, we artificially re- normalized 
all the profiles to equal residual intensities {r^ = 0.47 at the line center) and measured 
equivalent widths of the central narrow portion [vc < 0.6) of the resultant Ha profiles. In 
these conditions, the variation of the Ha core is attributed only to the rotationally modulated 
low-field component. As hoped, we find a close correlation between the EWcore measured in 
this way and the value of 15^1 measured by Maxted. This relation is shown in Fig. 4. 



^For distinctness we label all magnetic observables related to the low- field component with the subscript 
"G" , assuming its large scale (Global ) geometry. Observables related to the strong-field component will be 
labeled with the subscript "Spot" or "S"' . 



- 6 - 



Finally, using the \Bg\ — EW relationship derived from the high-resolution spectra as 
illustrated in Fig. 4, we inferred the field modulus \Bg\ associated with each of the low- 
resolution spectra (see Table [2]). 



4.2. The strong-field component 

In order to measure the magnetic field modulus \Bspot \ of the strong- field component we 
deblended the Ha profile by means of a simultaneous fit of five Gaussian profiles (three central 
profiles used to fit the Ha core, and two satellite gaussians to reproduce the strong-field 
Zeeman pattern). This method enabled us to reproduce the Zeeman splitting of the strong- 
field component and the corresponding magnetic field strength in those spectra where the 
strong-field spectral features are seen. The method also allows us to estimate the projected 
fractional area S of the strong-field area on the disc. Reconstructing by gaussians and 
extracting the Zeeman pattern of the strong field component from the observed Ha line 
profiles, we determined S", the fraction of the flux absorbed by the strong-field pattern 
relative to the total Ha absorption. The method is rather rough and can be considered as a 
first-guess approximation that is necessary for the analysis described below. A more realistic 
calculation of the size of the strong-field area is performed in Sect. 8 where we model the 
spectra. The results (|-Bspot| and S) are presented in Table [31 5* is given in per cent of the 
disc area. 



Mean longitudinal magnetic field of WD1953-011 



From the Stokes / and V spectra obtained with t he VLT and BTA we determined 
longitudinal fields through the weak- field approximation (jAngel et al. 1973 ) modified to the 
analysis of the two-component Stokes V spectra: 



m~(i-S)Bt(^)^^5i^+sv(AW (4) 

where Bq is the longitudinal field of the low- field component, describes the gradient of 
the fiux profile, S is the relative area of the spot projected on the disc, Aq is the Ha rest 
wavelength and V{\)spot is the Stokes V profile from the strong-field component observed 
in the Ha wings. 

In Eq.(4) we effectively separate the disc into two equivalent areas with different averaged 
magnetic field strengths. The first term in the equation describes the weak-field area and the 
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second term is attributed to the strong- field component. The first term is used in the usual 
manner according to which the flux and its gradient are taken directly from the observed 
spectra. (Inaccuracies due to the presence of the strong-field features in the wings are 
comparatively weak: these features are located quite far from the line core and do not affect 
the narrow central low-field polarization profile). However, circular polarization V{X)spot 
from the strong-field component (the second term in the equation) cannot be fitted in the 
same way. 

In order to fit V{X)spot and estimate the longitudinal field of the strong-field area, we 
compute: 

TA/\\ H^) Spot ~ H^) Spot /r-N 

V{X)spot = ^TV\ 



where the flux /(A) is the observed Ha flux profile, and I{\)spot ^"^^ H'^)spot 1^^" 
and right-hand polarized parts of the Ha profile from the strong-field equivalent area of the 
disc. In the observed polarization spectra I{X)spot I{X)spot mixed with fluxes from 
the weak-field equivalent area and therefore can not be extracted directly. However, we may 
estimate them with some simplifications. 

Individually, I{X)spot H^)spot Zeeman-split profiles of the circularly polarized 
satellite a components. Due to the fact that the a- component is absent in I{\)spoty 
the 0"+ component is absent in I{\)spoti their centers of gravity are displaced, indicating the 
presence of the longitudinal field from the strong-field area. Their difference provides the 
non-zero circular polarization (Eqn. 5). 

Because I{\)spot ^"^^ H^)spot ^'^^ resolved in the total left- and right-circularly po- 
larized observed fluxes, the determination of their true shapes requires detailed modeling the 
field geometry. However, as a first-guess approximation we may describe them by simulat- 
ing an equivalent "mean" Zeeman-broadened Ha profile magnetically displaced to the left- 
and right-sides from the rest wavelength [^. With this simplification only two parameters - 
Zeeman broadening and their magnetic displacement due to the averaged longitudinal field 
from the spot should be varied to reproduce the observed circular polarization. 



^The low resolving power of the FORSl and UAGS, as well as unresolved circular polarization features 
attributed to the strong-field area, enable us to consider the problem in terms of Zeeman broadening instead 
of detailed analysis of the strong-field Zeeman pattern. 
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In the simulation procedure we may use a template Ha profile, artificially broadened to 
unresolved Zeeman patterns typical for I{X)spot H^)spot magnetically displaced. 
This template can be taken from a zero magnetic field solution for the atmosphere of 
WD 1953-011 or from the observed spectra. For example, assuming the pressure-temperature 
conditions in the spot area to be similar to conditions in the other parts of the white dwarf's 
surface we may choose as the template one of the observed weak-field Ha profiles (obtained 
at those moments when the spot is not seen). In our analysis we proceed this way. 

Thus, to simulate I{X)spot H'^)spot order to fit the circular polarization (5) from 
the strong-field area and estimate its longitudinal magnetic field we used the following iter- 
ative method: 



• Step-1: We construct the reference weak-field "template" Ha profile from the observed 
I-profiles obtained at those rotational phases where the strong-field Zeeman pattern is 
not seen. 



Step-2: We artificially broaden the template profile by a gaussian filter with an arbi- 
trary half-width to an unresolved strong-field Zeeman pattern and displace the result 
by the magnetic displacement factor AA to the shorter / longer wavelengths to estimate 
the I{X)^spot and /(A)fp„, profiles. 



\R 



Step-3: Varying the magnetic broadening of the estimated profiles I{\)spot and I{^)spot 
and their Zeeman displacement we finally fit the strong-field circular polariz ation (5) in 
the w orking equation (4). The displacement found AA = 4.67- IO^^^B^Xq flLandstreet 



1980l ). gives an estimate of the longitudinal field B^. (In other words, taking 5* mea- 
sured from the Stokes I spectra and presented in Table [3] we simultaneously fit the 
observed combined Stokes V (4) varying Bq and circular polarization (5) of the strong- 
field component, where -B^ is one of the parameters.) 

In the fit procedure, the associated error bars ar e obtained using the Monte Carlo 



modeling method presented by lSchmidt fc SmithI (119941 ). An example of the fit is presented 



in Fig. O The results are collected in Table IH 

This method gives quite robust estimates of the longitudinal magnetic field of the low- 
field component. In the case of the strong- field component, the real intensities of the fields 
could be slightly over or underestimated due to the simplifications described above. For 
these reasons, estimates of the strong-field component given here could be considered to be 
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approximate. As an alternative, the tw o-component circula r polarization spectra could be 
analyzed by using Zeeman tomography (lEuchner et al.ll2002l . for instance). To provide more 
precise modeling, below (Sec. 8) we analyze our data again in the framework of simplified 
Zeeman tomography. 



6. Period determination 

To search for the star's rotation period we used the equivalent widths EWcore of the 
Ha core determined in Sect. 4 . This observable is the most sensitive indicator for the 
determination of the r otation period. To det ermine the rota tion period we a pplied the 
Lafler-Kinman method (ILafler fc Kinmanlll965l ). as modified by iGoranskyI (120041 ). Analysis 
of the power spectrum of the data revealed a signal indicating a probable period between 
1.4 and 1.5 days (Fig . [6]). This is consistent with t he period estimate (P ~ 1.45 days) given 
bv lWade et aP J2003h and ISrinkworth et all J2005h . 



Detailed study of the periodogram showed that the most significant sinusoidal signal 
corresponds to a period P = 1.4480 ± 0.0001 days. Other peaks are located around 1.447 
days and 1.442 days. An examination of these periods reveals distorted, non-sinusoidal 
signals and we do not consider these periods further. 

The phase variation of the Hct core equivalent widths EWcore derived with this period 
is presented in Fig. [71 The derived period shows a very good agreement among all the 
observations taken from different telescopes (the VLT, BTA and AAT) . For the minimum of 
EWcore we obtain the following ephemeris: 

JD = 2452048.801 ± 0.03 + 1^^4480 ± 0.0001 E 



The corresponding phase curves of the mean field modulus \Bg\ and longitudinal field 
\Bq\ of the weak-field components are presented in Fig. [HI The phase curves are almost 
sinusoidal, and symmetric about the values of about \Bg\ = -1-87 kG and \Bq\ = — 43kG. 
The modulus of the weak-field component varies from -|-77 ± 1.5 kG to -f97 ± 1.5 kG; the 
longitudinal magnetic field shows variation in the range —39 ± 2 kG to —47 ± 2 kG. 

As one can see in Fig. [8] the behavior of the longitudinal magnetic field and modulus of 
the weak-field component suggest the field geometry to be a simple, low-order poloidal field, 
which supports our view that the period is correct. Below we use this period in analyzing the 
magnetic field morphology of WD 1953-011. At the same time it is impo rtant to note that 
anoth er estimate of the rotational period (P = 1.4418 days) obtained by iBrinkworth et al. 



( 120051 ) from differential photometry of this MWD is similar to, but formally different from. 
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our result. Following the next section, where we establish the magnetic field morphology of 
the low- field component, we will discuss this difference in more detail. 



7. Modeling the weak-field component of the magnetic field of WD 1953-011 

To verify that the behavior of the weak- field component of WD 1953-011 is consis- 
tent with a nearly d ipolar geometry, we have followed the schematic method proposed by 



Landolfi et al.l (119971 ). This method has already been described and applied to establish the 



magne tic field morphology of the weak field white dwarf WD0009+501 (see IValyavin et al. 



( 120051 ) for details). For this reason here we do not explain all the modeling details, but 



restrict ourselves to the presentation of the results. 

In this paper we model the phase-resolved measurements of the mean longitudinal field 
and mean field modulus of the weak-field component within the framework of a pure dipole 
and dipole-|-quadrupole field. The phase-resolved observables for the weak-field component 
which we use as input data are obtained by binning measurements in phase and averaging. 
The binned data are presented in Table [51 

The dipole or dipole plus quadrupole models depend on the following 10 parameters: 

- Bd and B^, the dipole and quadrupole strength, respectively; 

- Ve, the stellar equatorial velocity; 

- i, the inclination of the stellar rotation axis to the line of sight; 

- f3, the angle between the dipolar axis and the rotation axis; 

- Pi and /?2, the analogues of P for the directions identified by the quadrupole; 

- 7i and 72, the azimuthal angles of the unit vectors of the quadrupole; 

- /o, the "reference" rotational phase of the model; 

- VeSini, the projected stellar rotation velocity. 

The angles i, P, Pi, P2 range from 0° to 180°, while 71, 72, /o range from 0° to 360°. 
The rotational period P = 1.448 and the limb-darkening constant u, which also affects the 
expressions of the magnetic observables, are taken as fixed. (Note that the pure dipole model 
would retain as free parameters only B^, Vg, i, P, and /q.) 
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For the stellar mass, iBragaglia et al.l (119951 ) gave the val ue of 0.844 Mp„ which together 
with known surface gravity of WD 1953-011 {log g = 8.412, IBragaglia et al.l (119951 )) corre- 
spond to a stellar radius of about 0.0095 Rq. This parameter and the period were then used 
to estimate the equatorial and projected velocities of the star. 

For t he limb-darken ir ig coe fficient, we adopted the value of u = 0.5. Note that, as dis- 
cussed by iBagnulo et al.l (120001 ) . the results of the modeling are only slightly influenced by 
the u value. The best-fit parameters are: 



A) Dipole 



i 


= 14° 


± 


10° 




= 14° 


± 


10° 


/o 


^ 352° 






Ba 


= 108 


± 


5kG 




= 0.33 


± 


0.05 km 


VeSini 


= 0.08 


± 


0.03 km s'^ 



B) Dipole + quadrupole 



i 


= 18° 


± 


10° 




8° 


± 


10° 


/o 


^ 357 






Pi 


= 22° 


± 


10° 




= 24° 


± 


10° 


71 


^ 77° 






72 


^ 243° 






5d 


= 178 


± 


30 kG 




= 233 


± 


30 kC 




= 0.33 


± 


0.05 km s 


VeSini 


= 0.1 


± 


0.05 km s 



Note that, as explained by iBagnulo et al.l (120001 ) . the available observations do not 
allow one to distinguish between two magnetic configurations symmetrical about the plane 
containing the rotation axis and the dipole axis. Such configurations are characterized by 
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the same values of B^, -Bq, fe, 7i, 72, fo, while the remaining angles are related by 

( h P, Pi, P2 ) 

( 180° - i, 180° - p, 180° - pi, 180° - /^s ) • 

Due to the fact that the spin axis angle is close to a pole-on orientation and due to 
the small number of available observables, the error bars on the derived quantities are fairly 
large. For the same reasons, there are some uncertainties in the determinations of all the 
parameters considered together. However, despite these weakness, we are able to obtain 
some conclusions about the most probable geometry of the white dwarf's global field. 

The best fit of the dipole+quadrupole model applied to the observations is shown by solid 
lines in Fig. O For comparison, the dashed line shows the fit obtained using the pure dipolar 
morphology. As is evident, the dipolar model does not reproduce the observations well. 
Examination of the reduced Xr statistics shows that the quality of the dipole+quadrupole 
fit {Xr = 0.22, 0.43 for longitudinal field and field modulus, respectively) is significantly 
better than the pure dipole fit {Xr = 1.9,2.2). We therefore conclude that the large-scale 
weak-field component of WD 1953-011 is better modeled by the superposition of a dipole 
and quadrupole components. 



8. High- field component 

8.1. Migrating magnetic flux tube? 

In contrast to the well-organized, nearly sinusoidal variation of the weak-field compo- 
nent, the phase behavior of the high- field structure exhibits a number of peculiar features 
that made it impossible to model these data as a simple low-order multipole: 

• According to the measurements of the Zeeman-split satellite spectral features in the 
Ha wings, the mean field modulus of the strong-field component does not show any 
noticeable variation during the star's rotation. (Due to rotation, we see variation of the 
fiux intensities from the strong-field area, but the corresponding Zeern an displacement 



is nearly constant.) The most likely explanation (iMaxted et al.ll2000l ) is that there is 
an area with a nearly uniformly distributed strong magnetic field. This explanation, if 
true, suggests that the strong-field component has a localized geometry and cannot be 
understood as a high-field term in the multipolar expansion of the star's general field. 
Averaging all the data we determine < |-Bspoi| >= 515 ± 7 kG. 
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The Zeeman pattern attributed to the strong-field component becomes visible at ro- 
tational phases = 0.25 — 0.7 and demonstrates variation in the flux intensities that 
suggests rotational variability of the projected effective size S of the magnetic spot. 
The projected area of the strong -fleld structure varie s from zero to about 12% of the 
disk, consistent with the study of iMaxted et al.l (120001 ). This observable can be used as 
an additional parameter to test the rotational period of the star. However, using this 
quantity to search for the period we did not flnd a regular signal at any period within 
the tested 5-year time base. Moreover, phasing the data with the magnetic ephemeris 
characterized by the rotational period of 1.448 days, the resultant phase curve of the 
spot size variation (Fig. 10) appears distorted in comparison to the well-organized be- 
havior of the weak-fleld component phased with the same ephemeris. We observe a 
small relative phase shift between the data obtained with different telescopes that may 
indicate a possible secular longitudinal drift of the strong-fleld component. 

The averaged longitudinal fleld of the strong fleld area varies from zero (when the 
area is invisible) to about 450 kG (Table Hj) which is comparable to the averaged 
mean fleld modulus (~ 515 kG) of the fleld. This suggests a deviation of the strong 
fleld component from any of low-order multipolar geometries for which the difference 
between the full vector and its longitudinal projection should much lar ger (for example, 
for a centered dipole fleld the difference should be at least 2.5 times, IStibbd (Il950l )). 



The last point suggests the presence of an essentially vertical orientation of the magnetic 
fleld lines relative to the star's surface, typical for local magnetic flux tubes in cool, convective 
stars (the Sun for example). If the geometry is a tube seen as a local magnetic spot in the 
photosphere, we may also expect the above-mentioned secular drift. To our knowledge 
and by an analogy to the Sun, such flelds are expected to show dynamical activity like 
migration over the star's surface and be associated with dark spots that might produce 
photometric varia bility of the star. Signiflcant photometric variability of WD 1953-011 has 
been established ( IWade et al.ll2003l : iBrinkworth et al.ll2005l ). However, in this paper we are 
unable to establish the association of the darkness and magnetic spots for reasons which we 
discuss below. 



8.2. Geometry and location of the magnetic spot 

Measuring the longitudinal magnetic fleld of the strong-fleld component as described in 
Sect. 5 we noted that the estimates of the longitudinal fleld of the strong-fleld area might 
be affected by our simplifying assumptions. In order to control these measurements we have 
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also directly modeled the observed polarization and intensity Ha spectra of the star without 
parametrizing the surface field components. In addition to the analysis of the Stokes / and 
V spectra, the observed Q and U spectra were also taken into consideration. As can be 
seen in Fig. 1, linear polarization Stokes Q,U signatures are detected only marginally at a 
few rotational phases. However, this information can also be used to constrain the magnetic 
geometry of the degenerate. Note that different magnetic geometries may produce similar 
Stokes V spectral features due to axial symmetry of circular polarization provided by the 
longitudinal projection of the field. Linear polarization restricts the strength and orientation 
of the transverse field that, together with circular polarization, makes it possible to resolve 
the geometry spatially. The observed Q and U spectra are mainly noise, but we may try to 
use them in terms of the upper limits. 

In our model we examined several low order multipolar magnetic field geometries, inte- 
grating over the surface elementary (taken at a single surface element) Stokes /, V, Q, and 
U spectra calculated for various field strengths and orientations of the magnetic field lines. 
The technique we used to calculate the elementary spectra deserves some special explanation. 



a) Simulation of the Stokes I,V,Q and U synthetic Ha spectra. 



Generally, accurate simulation of the split Balmer profiles in spectra of strong magnetic 
white dwarfs requires detailed computations of the main opacity sources under the influence 
of strong magnetic fields. To our knowledge, these computations have not yet been tabulated 
for practical use. For this reason, a self-consistent solution of the transfer equation for the 
line profiles in spectra of strong-magnetic white dwarfs cannot be performed without special 
consideration of additional parameters (related , for e xample, to the Stark broadening in 
the presence of a strong magnetic field, iJordaru (119921 )). However, in case of the weak- field 
degenerates we may restri c t our self to a zero- field so l ution similar to that presented by 
Wickramasinghe fc Martini (119791 ) or by lSchmidt et al.l (119921 ). The method assumes that if 
the Stark broadening dominates the line opacity, the total opacity can be calculated as the 
sum of individual Stark-broadened Zeeman components. The Stark broadening is suggested 
to be taken as "non-magnetic" in this case. Under these simplifying assumptions we may 
simulate the local, elementary Zeeman spectra using one of the following two ways: 



i to compute the transfer equation for all Stokes parameters at given strength and ori- 
entation of a local magnetic field line calculating the Ha opacities as described, or 
(alternatively) 



ii to select a "template" Ha profile typical for a zero-field white dwarf with the same 
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pressure-temperature conditions as in WD 1953-011 for construction of the elementary 
Zeeman spectra. (In other words, we may try to construct from this template profile 
individual Zeeman vr— and a— components, parametrizing their magnetic displace- 
ment and relative intensities, and additively combine them to obtain the elementary 
/, 1/, g, U Ha profiles.) 

The first, direct method of atmospheric calculations for WD 1953-011 requires special 
theoretical tools which are outside the scope of this observational paper. The second, simpli- 
fied method, which we will use, seems to be rather rough due to the fact that the fluxes from 
the individual vr— and a— components obtained by using the zero-field template profile are 
generally not additive (whereas their corresponding opacities can be added in the transfer 
equation). Nevertheless, in the linear guess approximation they can be taken as additive 
and the method can also be applied. Besides, testing this method on some standard, well- 
studied magnetic Ap/Bp stars we have obtained satisfactory results modeling their observed 
polarizations. This allowed us to conclude that the method is reasonably accurate. 

Thus, to calculate the elementary /, V, Q, U Ha profiles we adopted the use of the zero- 
field Ha template profile which was constructed from the observed I-spectra obtained at 
those moments, when the strong-field Zeeman pattern is not seen. The Stark parts of the 
profile were obtained by averaging the / profiles at the rotation phases and 0.91 (see 
Fig. 1) in which the strong- field features are not seen. The central "zero- field" Doppler 
profile was adopted to reproduce in the model procedure the observed low-field magnetic 
broadening of the Ha cores at phases and 0.91. The necessary individual profiles of 
the non-displaced vr— and displaced a— components were obtained by entering the normal 
Zeeman displacements according to the orientation of the local magnetic field: the circularly 
polarized a— components are displaced according to the longitudinal projection of the local 
magnetic field, and the linearly polarized a— components are displaced by the transverse 
field. 

In order to model the polarization Ha profiles, the relative intensities of the central vr— 
and displaced o"± components were computed from the prjection of t he loc a l field vector 



onto the plane on the sky an d on the line of s ight as described by lUnnd (119561 ) (for a 



qualitative explanation see also iLandstreetl (jl980l )). The final intensities of the vr— and a 



components were obtained by renormalisation such that the total sum of the fiuxes from all 
the components be equal to the fiux from the zero-field template Ha profile. 

Finally, the elementary vr— and a— components were combined to construct from them 
the elementary /, V, Q, U spectra by simulation of the ordinary and extraordinary beams 
given by a polarimetric analyzer. For example, simulating the Ha Stokes V profile, all the 
components except the circularly polarized a— components are equally distributed between 
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the beams. The circularly polarized 0"+ components are absent in one of the beams, and the 
oppositely polarized (T_ components are absent in the other beam. The final Stokes V Ha 
profile was obtained by subtraction of the ordinary from extraordinary beams and devision 
of the result by the total fiux. The Q and U spectra were obtained in a similar way: the 
linearly polarized a components and central tt component are distributed between the beams 
according to the projection of the local magnetic field onto the plane of the sky. 

After the determination of the elementary J, V, Q, U spectra given by a magnetic geom- 
etry in all surface elements, we finally integrated and averaged them over the disc. For the 
limb-darkening coefficient we adopted the value of u = 0.5. 



b) Modeling the field geometry in WD 1953-011 by simulation of the observed polarization 
spectra. 

In Sect. 5) we have concluded, that estimates of the longitudinal magnetic field based on 
the weak-field approximation are accurate enough to make it possible to model the geometry 
of the low-field component separately from the strong-field component in the manner as 
demostrated in Sect. 7. For this reason, and in order to reduce a number of variables, we use 
those results (obtained in Sect. 7) as input and non-changeable parameters in the tomography. 
We just note, that modeling the observed spectra obtained at those time moments, where 
the spot is not seen, we have confirmed, that the observations (Stokes-V spectra) can be 
better fit by the dipole+quadrupole geometry of the low-field component with parameters 
performed in Sect. 7 ( case B). Examination of the pure dipole model (case A) gives no 
satisfactory results and we do not use this case here. 

Modeling the strong-field area as an additional harmonic in the low-order (lower than 
octupole) multipolar expansion we were unable to reproduce the observations. The fit does 
not provide the necessary contrast in the observed Zeeman patterns at those phases where the 
polarization and intensity spectra demonstrate the weak- and strong-field Zeeman features 
together. However, assuming the strong-field component to be concentrated into a localized 
area having maximum projected size of about 12 per cent of the disk, the strong-field Zeeman 
spectral features can be well-reproduced with a n average magnetic field of 550 ± 50 kG. 



Practically the same result has been obtained by lMaxted et al.l (120001 ). 



To model the strong-field area we tested two simplified localized geometries: a "contrast 
spot" with a homogeneously distributed, essentially vertical magnetic field, and a "sagittal" 
geometry with a strong vertically-oriented central magnetic field, that smoothly decreased 
to zero at the spot edges. Generally, both geometries are able to describe the Stokes / and V 
spectra with more or less acceptable accuracy. The first model, however, does not provide us 
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with a good fit of the Q and U spectra due to the presence of the sharp (and non-physical) 
jump of the field intensity at the edges of the strong-field area. For this reason we do not 
discuss this case in detail. 

The "sagittal" geometry of the strong-field area was constructed by using a modified 
model of a centered dipole: about 45% of the spot's area (central parts) have the dipolar 
distribution with polar field Bp = +810 kG at the center. The remaining 55% of the external 
dipolar field is artificially modulated to have a gradual decrease to zero at the edges of the 
area. This model provides a good fit of the Stokes /, V spectra and reasonable reproduction 
of the linear polarization Q-,U-spcctra, as shown in Fig. 11 where we also illustrate the 
tomographic portrait of the white dwarf's magnetosphere. 

Despite the fact that we obtain such a good agreement of the "sagittal" geometry with 
the observables from the strong-field area, we do not claim that this geometry is fully correct 
in all details (for example, the model does not control conservation of magnetic fiux). Similar 
to the case of the weak-field component, the most natural way to study the strong-field area 
is to describe it as a strong-field feature resulting from the superposition of several high-order 
harmonics in the multipolar expansion. At this time we are unable to study this term using 
any combination of the first hamonics higher than octupole, but we do not exclude that the 
use of the highest terms of different polarities and intensities will resolve the problem. 

However, this result clearly demonstrates a qualitative difference in the morphologies of 
the strong-field area and the global field of the white dwarf. From the model we estabhsh with 
a very high probability that the strong-field area has a localized structure with essentially 
vertical orientation of the magnetic field lines. The physical size of the area is about 20% 
of the star's surface giving maximum 12% projection on the disc. The spot is located at 
an angle of about 67° with respect to the spin axis, providing a maximum longitudinal field 
strength of about 400 kG. These results are in good agreement with the measurements of the 
longitudinal field of the spot. 

9. Discussion 

We have presented new low-resolution spectropolarimetric observations of the magnetic 
white dwarf WD 1953-011. From these observations and observations of previous authors 
we have determined the star's rotation period, mean longitudinal field, mean field modulus, 
and surface field morphology. Let us finally summarize these results. 

1) Our present picture of WD 1953-011 consists of a MWD with relatively smooth, low- 
field global magnetic field component, and a high-field magnetic area. 
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2) The low-field component demonstrates regular periodicity with period P = 1.4480 ± 
0.0001 days. We interpret this as the rotational period of the white dwarf. The long- 
term stability of the surface and longitudinal magnetic fields of the low-field component 
enable us to interpret this component as a fossil poloidal magnetic field consisting of 
dipolar and quadrupolar harmonics with the following basic parameters: 

1) the inclination of stellar rotation axis i = 18° ± 10°; 

2) the angle between the dipolar axis and the rotation axis (3 = 8° ± 10°; 

3) the dipole strength = 178 ± 30 kG; 

4) the quadrupole strength Bq = 233 ± 30 kG. 



3) The strong-field component exhibits a peculiar localized structure. The mean field 
modulus of the spot \B,^p„f\ is estirnated t o be 515 ± 7 kG, which is consistent with the 



results presented by iMaxted et al.l ( l2000l ). The longitudinal magnetic field of the spot 
varies with rotational phase from < 300 kG to about 400 kG. Comparing the mean 
field modulus with the maximum longitudinal field we suggest that the geometry of the 
high-field spot may be similar to a magnetic flux tube with vertically-oriented magnetic 
field lines. The spot is located at an angle of ~ 67° with respect to the spin axis. 



Our results suggest that the magnetic field of WD 1953-011 consists of two physically 
different morphologies - the fossil poloidal field and an apparently induced magnetic spot. 
To our knowledge, fossil, slowly decaying global magnetic fields are organized in a nearly 
force-free poloidal configuration §. In contrast, if the suggested vertical orientation of the 
magnetic field lines in the spot is correct, the uncompensated "magnetic pressure"^ of such 
a localized field may dominate against the tension causing a strong impact on the pressure- 
temperature balance in the photosphere of the degenerate. This may produce a temperature 
difference between the strong-field area and other parts of the star's surface. As a result 
we may expect rotatio nall y-modulated photo r aetric variability of WD 1953-011 (as observed 
by I Wade et al.l (120031 ) and iBrinkworth et al.l (120051 ) ). For these reasons (and by analogy to 
sunspots) such fields might be unstable if not supported by other dynamical processes such 
as differential rotation, and may therefore exhibit secular drift with respect to the stellar 
rotation axis. 



■^According to the basic properties of the MaxweU stress tensor (see, for instance, iParkerj (|l979f )l the 
magnetic field B creates in the atmospheric plasma an isotropic pressure and tension directed along 
the magnetic lines of force. While neighboring lines of force of a magnetic field try to expand due to the 
magnetic pressure, tension tends to compensate for this effect. The force- free configuration is possible only 
when the gradient of the magnetic pressure is fully compensated by the tension forces. 
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In the above context we note that significant photometr ic variabihty of WD 1953-011 has 
been established ( jWade et al.ll2003l : iBrinkworth et al.ll2005l ). Also, remarkably, the authors 
had encountered problems analyzing the periodicity of the variable differential flux. When 
individual epochs of their photometric data are phased according to the rotation period of 
about 1.45 days, the resultant folded lightcurves are smooth and approximately sinusoidal. 
However, they had diffic ulty obtaining an a cceptable fit to all epochs of photometric data con- 
sidered simultaneously (IWade et al.ll2003l ). Besides, periodograms (IBrinkworth et al.ll2005l ) 
obtained separately for their 7 individual observing runs indicate a significant spread in the 
period distribution. The individual peaks are stochastically distributed around a rotation 
period P ^ 1^^.45 from P ^ 1^^.415 to P ^ 1*^.48, also suggesting a probable phase shift from 
epoch to epoch with characteristic times from tens to hundred of days. Combining all the 
data, they establish their version ot the rotation period P = 1. '^441769(8) days. This period 
is significantly different from the period derived by us from the behavior of the global field 
of the star (P = 1.448 ± 0.0001 days). 

Indirectly, these facts suggest a physical relationship between the darkness and magnetic 
spots, and their possible secular migration. Unfortunately, our spectroscopy and the available 
photometry were obtained at different epochs, making it impossible to study this relationship 
in this paper. Examination of these problems will be among the goals of our further study 
of this magnetic degenerate upon carrying out the necessary simultaneous photometric and 
spectral observations of this star. 
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Table 1: Spectral and spectropolarimetric observations of WD 1953-011: column 1 is the 
Julian Date of the midpoint of the observation, column 2 is the exposure time, and column 3 
reports the telescope used in the observa tions (AAT indicates the high resolution spec- 
troscopy presented by Maxted et al. ( 200ol )). For data obtained with the VLT the exposure 
times are presented for the three consecutive I iV / Q / U modes of observations (in this 
case the midpoint corresponds to observations of the Stokes /, V parameters). 



JD 


Exp (sec) 


Stokes 


Telescope 


2450676.955 


600 


/ 


AAT 


2451391.948 


600 


/ 


AAT 


2451391.955 


600 


/ 


AAT 


2451391.962 


600 


/ 


AAT 


2451392.059 


1800 


/ 


AAT 


2451392.957 


1800 


/ 


AAT 


2451393.066 


1800 


/ 


AAT 


2451393.106 


1800 


/ 


AAT 


2451393.943 


1200 


/ 


AAT 


2451393.958 


1200 


/ 


AAT 


2451393.973 


1200 


/ 


AAT 


2451393.988 


1200 


/ 


AAT 


2451394.003 


1200 


/ 


AAT 


2452048.801 


840/1200/1200 


I,V/Q/U 


VLT 


2452048.893 


840/1200/1200 


I,V/Q/U 


VLT 


2452076.671 


840/1200/1200 


/, V/Q/U 


VLT 


2452076.883 


840/1200/1200 


I,V/Q/U 


VLT 


2452078.722 


840/1200/1200 


I, V/Q/U 


VLT 


2452078.879 


840/1200/1200 


I, V/Q/U 


VLT 


2452079.672 


840/1200/1200 


I, V/Q/U 


VLT 


2452079.892 


840/1200/1200 


I, V/Q/U 


VLT 


2452087.621 


840/1200/1200 


I, V/Q/U 


VLT 


2452087.670 


840/1200/1200 


I, V/Q/U 


VLT 


2452087.722 


840/1200/1200 


I, V/Q/U 


VLT 


2452087.768 


840/1200/1200 


I, V/Q/U 


VLT 


2452505.290 


3600 


i,v 


BTA 


2452505.327 


3600 


i,v 


BTA 


2452505.360 


3600 


i,v 


BTA 


2452505.397 


3600 


i,v 


BTA 
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Table 2: Determinations of the mean modulus \Bg\ of the weak-field component. Column 1 
is the Julian Date, Col. 2 and Col. 3 are the equivalent widths of the Ha core (EWcore) and 
associated error bar, Col. 4 and Col. 5 are the inferred field strength \Bg\ and its error bar 
a(kG), Col. 6 is the telescope used. Uncertainties at the measured equivalent widths are 
calculated as a noise fraction of the flux (due to Poisson noise) in the total flux under the 
line profile. The mean field modulus and its u ncertainty obt a ined f rom the high-resolution 
spectroscopy with the A AT are taken from iMaxted et al.l (|2000| ). Uncertainties at the 
calibrated field strengths (observations with the VLT and BTA) result from regression 
errors in the EWcore^lBcl relationship shown in Fig. 4. 



JD 


EWcore 


a 


\Bg\ (kG) 


a{kG) 


Telescope 


2450676.955 


1.040 


0.008 


91 


5 


AAT 


2451391.948 


1.073 


0.008 


93 


4 


AAT 


2451391.955 


1.105 


0.016 


100 


4 


AAT 


2451391.962 


1.045 


0.008 


93 


4 


AAT 


2451392.059 


1.068 


0.012 


93 


2 


AAT 


2451392.957 


0.905 


0.012 


83 


1 


AAT 


2451393.066 


0.943 


0.012 


80 


2 


AAT 


2451393.106 


0.933 


0.012 


83 


1 


AAT 


2451393.947 


1.008 


0.020 


92 


3 


AAT 


2451393.958 


1.013 


0.016 


87 


2 


AAT 


2451393.973 


0.945 


0.016 


84 


2 


AAT 


2451393.988 


0.935 


0.020 


84 


2 


AAT 


2451394.003 


0.981 


0.020 


83 


2 


AAT 


2452048.801 


0.901 


0.012 


80 


3 


VLT 


2452048.893 


0.917 


0.012 


81 


3 


VLT 


2452076.671 


1.012 


0.012 


89 


3 


VLT 


2452076.883 


1.057 


0.016 


93 


3 


VLT 


2452078.722 


1.013 


0.016 


89 


3 


VLT 


2452078.879 


0.939 


0.012 


83 


3 


VLT 


2452079.672 


1.052 


0.016 


93 


3 


VLT 


2452079.892 


1.102 


0.012 


97 


3 


VLT 


2452087.621 


0.952 


0.016 


82 


3 


VLT 


2452087.670 


0.939 


0.016 


83 


3 


VLT 


2452087.722 


0.911 


0.016 


78 


3 


VLT 


2452087.768 


0.921 


0.016 


79 


3 


VLT 


2452505.290 


0.978 


0.020 


84 


4 


BTA 


2452505.327 


1.000 


0.020 


88 


4 


BTA 


2452505.360 


1.028 


0.020 


91 


4 


BTA 


2452505.397 


1.034 


0.020 


91 


4 


BTA 
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Table 3: Determinations of the mean modulus |-Bspot| of the strong-field component. Column 1 
is the Juhan Date, Col. 2 and Col. 3 are relative area of the spot on the disc S (in percent 
of the full disc area) and associated error bar obtained as a noise fraction of the flux in the 
total flux under the strong-field sateUite features. Col. 4 and Col. 5 are the magnetic field 
strength |Sspot| and its error bar a (kC) obtained as uncertainty in the determination of the 
satellite positions deblended by Gaussians. Col. 6 is the telescope used.) 



JD 


SiVo) 




spot \ \'^^ J 


a{kG) 


Telescope 


2450676 955 


13.1 


0.6 


521 


40 


AAT 


2451391 948 


9.5 


0.6 


513 


30 


AAT 


2451391 955 


12.0 


1.2 


495 


30 


AAT 


2451391 962 


11.7 


0.6 


494 


30 


AAT 


2451392 059 


12.0 


0.6 


527 


15 


AAT 


2451392 957 


0.9 


0.9 


iTivi si Hip 




AAT 


2451393 066 


0.6 


0.9 


in VI si nip 




AAT 


2451393 106 


0.7 


0.9 


invi^iHlp 




AAT 


2451393 947 


5 4 


1 5 


invisi nip 

±11 V loi iL'iv;; 




AAT 


2451393 958 


2.4 


1.2 


invisiVilp 




AAT 


2451393 973 


0.6 


1.2 


invi^iVilp 




AAT 


2451393 988 


0.7 


1.5 


in visiVilp 

111 V lOl > J i- V. 




AAT 


2451394 003 


0.7 


1.5 


invi si Hip 

111 V lOl k-Zlv/ 




AAT 


2452048.801 


0.6 


0.9 


invisible 




VLT 


2452048.893 


0.6 


0.9 


invisible 




VLT 


2452076.883 


7.8 


1.2 


520 


15 


VLT 


2452078.622 


3.2 


1.2 


invisible 




VLT 


2452078.879 


0.5 


1.2 


invisible 




VLT 


2452079.672 


8.4 


1.2 


529 


30 


VLT 


2452079.892 


12.3 


0.9 


511 


15 


VLT 


2452087.621 


0.6 


1.2 


invisible 




VLT 


2452087.670 


0.6 


1.2 


invisible 




VLT 


2452087.722 


0.7 


1.2 


invisible 




VLT 


2452087.768 


0.6 


1.2 


invisible 




VLT 


2452505.290 


10.4 


1.2 


500 


35 


BTA 


2452505.327 


12.3 


1.2 


492 


40 


BTA 


2452505.360 


12.3 


1.2 


524 


35 


BTA 


2452505.397 


10.8 


1.5 


502 


45 


BTA 
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Table 4: Determinations of the longitudinal magnetic field of the low- and high-field compo- 
nents of WD 1953-011. Column 1 is the Julian Date, Col. 2 and Col. 3 are the longitudinal 
field of the low-field component and associated error bar, Col. 4 and Col. 5 are the deduced 
longitudinal magnetic field of the strong- field component and its error bar ("no" means 

"l)c4()w the (letectioii levc4"). Col. 6 is the telesc()])e used.) 



JD 




a (/,-G) 


{kG) 


a {kG} 


OBS 


2452048.801 


-41.5 


1.5 


no 




VLT 


2452048.893 


-39.6 


1.6 


no 




VLT 


2452076.883 


-41.0 


1.6 


430 


70 


VLT 


2452078.722 


-42.9 


1.8 






VLT 


2452078.879 


-42.2 


1.7 


no 




VLT 


2452079.672 


-41.9 


1.6 


360 


60 


VLT 


2452079.892 


-46.8 


1.7 


460 


60 


VLT 


2452087.621 


-41.5 


1.6 


no 




VLT 


2452087.670 


-39.8 


1.7 


no 




VLT 


2452087.722 


-40.1 


1.7 


no 




VLT 


2452087.768 


-40.1 


1.5 


no 




VLT 


2452505.290 


-46.2 


2.3 


440 


80 


BTA 


2452505.327 


-44.8 


2.3 


450 


80 


ETA 


2452505.360 


-45.0 


2.5 


no 




BTA 


2452505.397 


-42.0 


2.7 


no 




BTA 
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Table 5: Phase- resolved observable magnetic quantities of the weak- field component of 
WD 1953-011. The first column is rotation phase obtained according to the magnetic 
ephemeris characterized by the rotational period of 1.448 days found here; the second and 
third columns are the mean field modulus or "surface magnetic field" \Bg\ and associated 
error bar; the fourth and fifth columns are the longitudinal field and its error bar. 








a 




a 




(kG) 


(kG) 


(kG) 


(kG) 


0.05 


81.5 


4 


-40.5 


1.3 


0.15 


82 


2 






0.25 


89 


3 






0.35 


94 


3 


-44 


1.5 


0.45 


93 


2 






0.55 


97 


3 


-47 


2 


0.65 


89 


3 






0.75 


86 


4 


-42.2 


2 


0.85 


83 


1 


-41 


1 


0.95 


78 


2 


-40 


1 
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Fig. 1.— Stokes IQUV Ra profile timeseries of WD1953-011, obtained using the FORSl 
spectropolarimeter at the ESO VLT. Phases correspond to the magnetic ephemeris obtained 

in this paper and arc expressed in part per mil at right. The thin lines represent the 
observations obtained at phase 0, and are reproduced to emphasise the variability of the 
Stokes profiles. 
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Fig. 2. — Stokes V at the Ha line obtained with the VLT (the solid line) and the BTA (the 
dashed line) . The spectra correspond to phases of the maximum visible strong-field Zeeman 
satellite features at the Ha wings {(f) ^ 0.3 in observations with the BTA and ^ ^ 0.5 in 
observations with the VLT: the observed phase shift is discussed in Sect. 8 and Sect. 9 of this 
study.) 
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Fig. 3. — The Ha profiles obtained at the VLT, AAT and BTA. High resolution spectra are 
convolved to the spectral resolution of the FORSEl and UAGS. The solid lines illustrate 
profiles at two extreme rotation phases at which the spot component is most clearly seen 
(the shallowest profile) and where the spot component is absent (the deepest profile). 
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Fig. 4. — The relationship EWcore ^ \Bg\ o btained from the con volved high- resolution 
spectra for which \Bg\ values are estimated by iMaxted et al.l (120001 ). The dotted hne is a 
linear fit of the relationship. 
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Fig. 5. — An example of the model technique as described in sect. 5. The thick solid line is 
the Stokes V observed spectrum containing strong circular polarization from the strong-field 
component; the dotted and dashed hues are modeled Stokes V spectra of the weak- and 
strong-field components respectively; the thin solid line is their sum. 
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Fig. 6. — Power spectrum of the magnetic field variations in WD 1953-011. 
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Fig. 7. — The phase curve of the equivalent widths at the Ha core phased with the 1.4480-day 
period. 




Fig. 8. — The magnetic phase curves of WD 1953-011 and their fits with the 1.448-day 
period. The upper plot illustrates variation of the field modulus of the weak-field component; 
the lower plot is variation of the its longitudinal magnetic field. 
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Fig. 9. — Observations and modeling of mean longitudinal field (bottom panel) and mean 
field modulus, or surface field (top panel). The dashed lines show a fit obtained with dipole 
model. The solid lines show the best-fit obtained by means of a dipole + quadrupole model. 



-33- 



15.0 
12.5 
10.0 
^ 7.5 
5.0 
2.5 
0.0 



0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 

Rotational Phase 

Fig. 10. — Phase-resolved projection S of the strong-field area: open circles illustrate the 
data from the AAT, filled circles and triangles are the observations with the VLT and BTA 
respectively. All the data have been phased according to the magnetic ephemeris obtained 
in Sect. 6. 
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Fig. 11. — Model fits (red lines) of the observed (black lines) Stokes IV QU spectra obtained 
in observations with the VLT. The corresponding tomographic portraits of the star's mag- 
netosphere and rotational phase (marked as Phase ) are presented at left. The strong- field 
area is shown by white. The magnetic field line of force are red lines. 
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Fig. 12. — continued. 



